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Abstract 



We consider the combined effects of a power law Levy step distribution 
characterized by the step index / and a power law waiting time distribution 
characterized by the time index g on the long time behavior of a random 
walker. The main point of our analysis is a formulation in terms of coupled 
Langevin equations which allows in a natural way for the inclusion of external 
force fields. In the anomalous case for / < 2 and g < 1 the dynamic exponent 
z locks onto the ratio f/g. Drawing on recent results on Levy flights in the 
presence of a random force field we also find that this result is independent of 
the presence of weak quenched disorder. For d below the critical dimension 
d c = 2f — 2 the disorder is relevant, corresponding to a non trivial fixed point 
for the force correlation function. 
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I. INTRODUCTION 



Anomalous diffusion arising from non trivial waiting time distributions, so-called contin- 
uous time random walk (CTRW), has been used to model a variety of physical phenomena 
0. For example, the dynamics of carrier diffusion and recombination in disordered media 
has been described in terms of a CTRW P-p]. Anomalous diffusion in an intermittent dy- 



namical system ||,[| and in a linear array of convection cells []T0| , [ll|1 has also been analyzed 
in terms of CTRW. 

Anomalous diffusion is also associated with power law step size distributions, the so-called 
Levy distributions The Levy flights (LF) generated by the Levy step distribution 



have been used to model a variety of physical processes such as self diffusion in micelle 
systems [|ll|] and transport in heterogeneous rocks 

In the present paper we consider the combined effect of a power law waiting time distri- 
bution and a power law Levy step size distribution, i.e., the case of continuous time Levy 
flights (CTLF). Since Levy flights lead to superdiffusive behaviour, whereas a power law 
waiting time distribution entails subdiffusive behaviour, the combination of the two yields 
an interesting description of anomalous diffusion in the general case 0. 

The key issue in the paper is the discussion of CTLF in terms of coupled Langevin 
equations for the position and the time. This approach is achieved by introducing as an 
intermediate variable the path or arc length along a particular trajectory. Besides the 
clearer physical interpretation afforded by Langevin equations as stochastic equations of 
motion such an approach also allows in a natural way for the inclusion of external fields, for 
example a drift force field. 

In section II we discuss the Langevin equation associated with power law step and waiting 
time distributions. Section III is devoted to the discussion of the associated Fokker-Planck 
equations. In section IV we discuss the probability distribution for the position of the 
walker as a function of time. In section V we carry out a scaling analysis in the absence of 
an external force field. In section VI we give a discussion and comment briefly on the role 
of an external quenched force field, drawing on recent results. 



II. LANGEVIN EQUATIONS 

Let us denote the step size distribution by ^(77), where r\ is the elementary microscopic 
step. We assume that Tr(r]) is normalized, i.e., / TT(r})d d i] = 1. Likewise, we denote the 
waiting time distribution by w(t), where w(t) is the probability of the walker of waiting the 
time interval r at a given position before performing the next step. Due to causality r > 
and the normalization of w{t) reads w(t)cLt = 1. 

Parametrizing the random walk in the continuum limit by means of the path param- 
eter or arc length s along the trajectory we have for the position of the walker, r(s), after s 
'steps', 

r( S )= f\(s')ds', (1) 

•J 

or the Langevin equation, 
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In Fig. 1 we have shown the parametrization of a particular trajectory. Equation (2) 
is readily generalized in the presence of an external drift force field F depending on the 
position of the walker. We obtain in this case the Langevin equation, 

-=F(r) + rj( S ). (3) 

Similarly, the total elapsed time after s 'steps' in the continuum limit is 

t{s)= [ S r(s')ds', (4) 
Jo 

implying the Langevin equation, 

The coupled Langevin equations (2), (3), and (5) constitute the present formulation of con- 
tinuous time random walk. We are here considering the separable case where the step distri- 



bution ir(rj) and the waiting time distribution w(r) are statistically independent [3,|i^,|i~8]| . 

In the case of a sharp waiting time distribution w(r) = 5(r — r ), corresponding to a 
fixed hopping rate l/r , Eq. (5) becomes deterministic and can be solved for s, i.e., t = r s, 
and we obtain by insertion the usual Langevin equation, 

^ = -F(r) + -n(t). (6) 

at t T 

In the general case of a non trivial waiting time distribution we must, however, discuss the 
coupled equations. 

We shall now focus on power law distributions for the waiting times and steps. For the 
step 77 we assume an instantly (in terms of s) correlated power law distribution, 

Tc(r])d d r] oc rf x ~ s dr}. (7) 

We assume an isotropic form characterized by the step index /. In order to ensure nor- 
malizability we have introduced a lower cut rj ~ a of the order of a microscopic length a 
and chosen / > 0. For / > 2 the second moment, (if) = J Ti{r])rf'd d rj ) is finite and a 



characteristic step size is given by the root mean square deviation J (r] 2 ). For 1 < / < 2 the 
second moment diverges but the mean step, (rj), is finite. In the interval < / < 1 the first 



moment diverges and even a mean step size is not defined ||12|| . In a similar way we assume 
for the waiting time r the power law distribution, 

w(T)dT oc T~ x ~ 9 dr. (8) 

Owing to causality r > 0. In order to guarantee normalizability we introduce a short time 
cutoff of order a microscopic time scale and choose g > 0. The distribution is characterized 
by the time index g. For g > 1 the first moment, (r) = / w(r)rdr, is finite, setting a well 
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defined hopping rate 1/(t). For < g < 1 the mean value (r) diverges and we cannot define 
a characteristic rate or time scale. 

In order to eventually discuss the anomalous diffusive characteristics of a random walk 
driven by the power law noises rj and r we must determine the distribution function P(r, t) 
and, in particular, the mean square displacement 

<r 2 (t)> ex Dt 2 '\ (9) 

where D is the diffusion coefficient and z the dynamic exponent. Since the random walk 
takes place in physical time t we are thus faced with the issue of eliminating the auxiliary 
path variable s labelling the walk. Owing to the stochastic nature of the Langevin equation 
(5) for t it is not possible to invert it and solve it for s and we shall instead turn to the 
associated deterministic Fokker-Planck equations. 



III. FOKKER-PLANCK EQUATIONS 

Limiting first the discussion to the force free case, i.e., F = 0, the 
probability distributions Px(t,s) and P 2 (t,s) associated with the Langevin equa- 
tions (2) and (5) are easily inferred fl9U2"0"| . From the definition P\(r, s) = 
(S(r — r(s))) and P2(t,s) = (5{t — t(s))), the Eqs. (1) and (4), and averaging according 
to the power law noises in Eqs. (7) and (8), we deduce the scaling forms, || 

^■•) = /w^' = ^" G >(^)' (10) 
PA S) = I wf'"'~'"" r ° = S_1/ " G2 (^) ■ 

By choosing a suitable scale for t and s we have fixed the coefficients of fc M and (— %uj) v to 
be unity. 

The scaling exponents \i and v depend on the step index / and time index g, respectively, 
characterizing the power law step size and waiting time distributions in Eqs. (7) and (8). 
For / > 2, i.e., the case of a finite mean square step, \i locks onto the value 2 and the scaling 
function G\ takes the Gaussian form characteristic of ordinary Brownian walk, G\[x) = 
exp (—x 2 ). This is a consequence of the central limit theorem which here leads to universal 
behaviour and defines the universality class of Brownian motion. For / < 2 the scaling 
exponent /x = /, and the scaling function G\ can only be given explicitly in terms of known 
functions for /i = 1 and \i — 1/2 (the Cauchy and Smirnov distributions, respectively, lfL7[l ). 
It is, however, easy to show that G\ — > const, for x — > and G\ — > for x — > oo. From the 
distribution in Eq. (10) we infer the scaling form for the mean square displacement of the 
walker in terms of the path variable s, 

(r\ s )) = J p^r, s)r 2 d d r oc s 2/ '\ (12) 

A similar discussion applies to P2(t,s). For g > 1, i.e., the case of a finite first moment, 
corresponding to a well defined hopping rate, the scaling exponent v locks onto 1 and the 



4 



scaling function G 2 (x) = 5(1 — x). This is again a consequence of the central limit theorem 
which leads to a universal time behaviour, i.e., a hopping rate. For g < 1 the scaling 
exponent v locks onto g, i.e., v = g and we obtain a non trivial scaling behaviour. For the 
mean square displacement in terms of the path variable we have, 

(t 2 (s)} = J P 2 (t, s)t 2 dt oc s 2/u . (13) 

In Fig. 2 we have depicted the scaling exponents /i and v as functions of the step index / 
and the waiting time index g, respectively. 

A simple heuristic argument || using Eq. (13) to infer that t scales like s 1 ^ and elimi- 
nating s in Eq. (12) yields the scaling relation, 

(r 2 (t)) oc t 2u/ ^ (14) 

and according to Eq. (9) the dynamic exponent, 

, = (15) 

As we shall discuss in section V this is indeed the correct result following from the scaling 
law for P(r, t), 

P(r,t)=t~fG(r/t^). (16) 

We conclude this section by writing down the Fokker-Planck equations following from Eqs. 
(10) and (11). Introducing the 'fractional' non local differential operators, 

rl d k 

e* kr F, (17) 



(27T) 

D" = I ^e'^i-iuy, (18) 

Z7T 



reflecting the long range Levy steps and waiting times, we have 

dPAr, s 



ds 

dP 2 (t,s) 
ds 



= VPi(r,s) (19) 
= D u P 2 (t,s). (20) 



In the Brownian case /i = 2 and v — 1 and V M reduces to the usual Laplace operator 
A describing ordinary diffusion, whereas D v becomes the first order differential operator 
-d/dt. 

In the presence of a force field F(r) we have, correspondingly, [|19[ 



dP f s ' S) = -V(F(r)P 1 (r, 8 )) + VPi(r, s). (21) 

Here the first term on the right hand side of Eq. (21) is the usual drift term due to the 
motion of the walker in the force field. 



5 



IV. THE DISTRIBUTION P(r,t) 



In order to calculate the probability distribution for the walker as a function of the 
physical time t we must eliminate the path variable s. In other words, we have to derive the 
distribution P 3 (s,t) since it then follows that P(r, t) is given by the relationship, 

P{r,t) = / dsP 1 {r,s)P 3 {s,t), (22) 

Jo 

due to the fact that the probability of the walker arriving at r in time t equals the probability 
of being at s on the path at time t multiplied by the probability of being at position r for 
this path length s, summed over all path lengths. 

In order to derive P 3 (s, t) we use the general expression for P(r, t) for arbitrary step and 



waiting time distributions given by Montroll and Weiss ||21|| . In Fourier space we have 



P(k,^) = (-^)- 1 1 (23) 
1 — w{Ld)7r{k) 

where 7r(k) and w{uj) are the Fourier transforms of the step and waiting time distributions, 
i.e., 

7r(k) = / e- ik - v 7c{ri)d d r] (24) 



w{u) = / e iuJT w(r)dT. (25) 



We notice that since w(r) = for r < due to causality Eq. (25) reduces to a Laplace 
transform and u is defined along a contour parallel with the x-axis in the upper half complex 
uj plane. 

For the power law distributions in Eqs. (7) and (8) we obtain in particular 7r(k) = 1 — 
for small k and w(u) = 1 — (—iuj) u for small u, i.e., to leading order in uo and k, the 
distribution 



P(k,u) = / . / , . (26) 

Inserting Eq. (10) for Pi(r, s) in Eq. (22) and requiring that P(r, t) is the Fourier transform 
of Eq. (26) it is easy to demonstrate that P 3 (s,t) is given by [0 

P 3 (s,t)= I ^-e-^H-iujy-h-^" 8 . (27) 

We have not found a simple physical argument leading to Eq. (27), i.e., the 'inversion' of 
Eq. (11), but notice that dP 3 /dt = dP 2 /ds. 
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V. SCALING ANALYSIS 



In the absence of the force field, i.e., for F = 0, it is an easy task to carry out a scaling 
analysis. The results are most easily deduced from Eq. (26). We have 

p (r f] _ [ duJ ^ .-Wkr {-M V ~ X 

P(r ' t] ~ J 2^ (27p e (-iu,)» + k»' (28) 

leading to the scaling form in Eq. (16) with scaling function 

_ d r duo d d k e -iu+ik 
G W- X J 2^ (2ttY (- tU j)[l + x-^/(-iujyy " ' 

The mean square displacement is given by 

(r 2 (t)) = J d d rP(r,t)r 2 oc t 2/z (30) 

and we find the dynamic exponent z — fi/v in agreement with the heuristic result in Eq. 
(15). 

Along the line fx = 2v we have z = 2 and CTLF have the same scaling characteristics as 
ordinary Brownian motion. The superdiffusive behaviour induced by the long range Levy 
steps is precisely balanced by the long waiting times. We notice, however, that the scaling 
function G depends on the scaling index fx = 2v\ only in the case fx — 2 and v = 1 do we 
obtain the Gaussian distribution. 

For fx > 2v we have z > 2, the Levy flights prevail, and we obtain superdiffusive be- 
haviour; correspondingly, fx < 2v implies z < 2, the long waiting times dominate, and we 
have subdiffusive behaviour. 

In Fig. 3 we have shown the different universality classes for CTLF. For z = 2 (fi = 2v) 
we have the universality class of ordinary Brownian motion. For fi > 2v we obtain the 
universality class (or classes) of anomalous super diffusion with an exponent z depending 
continuously on the ratio of the microscopic exponents / and g. Similarly, for fi < 2v we 
have the universality class of anomalous subdiffusion. 



VI. DISCUSSION AND CONCLUSION 

In the present paper we have discussed the combined effects of an algebraic waiting time 
distribution and an algebraic Levy type step distribution on the motion of a random walker. 
In order to include external force fields we have formulated this analysis in terms of a set of 
coupled Langevin equations. In the absence of force fields a simple scaling analysis shows 
that the dynamic exponent z characterizing the long time behaviour of the mean square 
displacement is given by the ratio z = fi/u, where the scaling exponents fx and v are related 
to the microscipic step and waiting time exponents {p, — f for / < 2 and v = g for g < 1). 
This dependence defines three universality classes: 1) Normal diffusive behaviour for ft = 2v, 
2) Anomalous superdiffusion for fi > 2u, 3) Anomalous subdiffusion for fi < 2v. 

In the presence of a quenched (time independent) Gaussian random force field and a 



well defined hopping rate, i.e., for v — 1, we have recently shown JT9| that z locks onto 
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fj, for / < 2. It is now easy to generalize this result to the case of a non trivial waiting 
time distribution with v < 1. Since the force field is independent of t we can directly apply 
the discussion in ref. [|19j to Eq. (3) treating the path variable s as the effective time. We 
conclude that in this case z is also given by fi/u for weak quenched disorder. The analysis 



in ref. ||19|| , leading to the critical dimension d c = 2/i — 2 below which a non trivial force 
correlation fixed point emerges can also be carried over to the present case. From a physical 
point of view it is clear that a nontrivial waiting time distribution has no effect since the 
quenched force field acting at position r simply 'waits' till the walker arrives. However, in 
the case of a time dependent random force field the waiting time distribution will interfere 
with the temporal force correlations and we have to treat the coupled Langevin equations in 
order to eliminate the intermediate path variable s. This interesting case will be considered 
in a forthcoming publication p2|. 
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FIGURES 




FIG. 1. Plot of a particular random walk parametrized by the arc length s 
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FIG. 2. In a) we show \i as a function of /. For / > 2 the exponent /U = 2 and we have ordinary 
Brownian walk; for / < 2 the exponent fi = f and we have Levy flights, leading to anomalous 
superdiffusion. In b) we depict v as a function of g. For g > 1 the exponent v = 1 and we have 
a well defined hopping rate for the walker; for g < 1 we have v = g and we obtain a subdiffusive 
behaviour. 
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FIG. 3. Plot of jJL versus v. Along the solid line jjl = 1v we have the universality class of 
ordinary Brownian motion with z = 2. In the region jjl > 2v we have the universality class of 
anomalous superdiffusion; in the region [i < 2v the class of anomalous subdiffusion. 
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